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ABSTRACT 


The  notion  of  a  generalized  finite  element  method  is  intro¬ 
duced.  This  class  of  methods  is  analyzed  and  their  relation  to 
mixed  methods  is  discussed.  The  class  of  generalized  finite  ele¬ 
ment  methods  offers  a  wide  variety  of  computational  procedures 
from  which  particular  procedures  can  be  selected  for  particular 
problems.  A  particular  generalized  finite  element  method  which 
is  very  effective  for  problems  with  rough  coefficients  is  dis¬ 
cussed  in  detail. 
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1 .  Introduction 

Although  the  general  ideas  involved  in  finite  element 
methods  (FEMs)  are  very  old  and  are  in  fact  closely  related  to 
classical  variational  principles  (see  [1  ],  e.g.,  for  a  history  of 
the  finite  element  method) ,  the  full  power  of  these  methods  is 
only  seen  in  their  implementation  as  computer  codes.  (For  a  list 
of  over  650  codes  for  finite  element  analysis  of  structural  me¬ 
chanics  problems  see  [  2  ]).  The  success  of  finite  element  methods 
is  due  to  several  factors. 

a)  They  allow  for  effective  implementational  procedures  based 
on  a  collection  process  which  leads  to  the  construction  of 
a  sparse  discretized  system  (stiffness  matrix)  . 

b)  Finite  element  methods  are  robust.  In  many  contexts  this 
robustness  is  a  consequence  of  underlying  variational 
principles  in  mechanics. 

c)  These  methods  allow  for  the  possibility  of  effective  post¬ 
processing^  as,  e.g.,  in  the  calculation  of  stresses  at 
arbitrary  points. 

d)  Many  finite  element  methods  have  a  physical  interpreta¬ 
tion  which  permits  an  intuitive  understanding  of  the  com¬ 
putational  process. 

From  the  point  of  view  of  present  day  computational  studies, 
properties  a),  b) ,  and  c)  (especially  a)  and  b) )  are  essential. 
Although  d)  is  important,  its  importance  has  diminished  because 
of  progress  in  the  mathematical  understanding  of  finite  element 
methods,  the  expanding  field  of  applications  for  these  methods. 
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and  the  progress  in  computer  technology. 

Since  FEM  s  are  increasingly  used  by  unsophisticated  users 
in  a  widening  variety  of  fields,  it  is  important  to  generalize 
the  notion  of  the  finite  element  method  as  far  as  possible  while  re¬ 
taining  the  basic  implementational  architecture  of  the  standard 
finite  element  methods  (see  a)  above),  so  that  a  particular  version 
can  be  selected  for  a  class  of  problems  so  as  to  be  maximally 
robust  and  to  allow,  if  possible,  for  effective  postprocessing. 

This  paper  addresses  the  problem  of  generalizing  the  notion 
of  the  finite  element  method  in  the  context  of  a  simple  one  di¬ 
mensional  problem.  We  restrict  ourselves  to  this  model  problem 
in  order  to  clearly  illustrate  the  ideas;  the  extension  to  two 
dimensional  problems  will  be  presented  in  a  forthcoming  paper. 

We  thus  consider  the  boundary  value  problem 

/r-(a(x)u')'  +  b(x)u  =  f , 0  <  x  <  1  ('  = 

Vu(0)  =  u(l)  =  0. 

Robustness  of  a  FEM  for  this  problem  will  mean  that  the 
method  performs  well  independently  of  the  smoothness  of  the  coef¬ 
ficients  a  and  b  and  the  source  term  f.  More  precisely,  we 
will  say  a  method  is  robust  with  respect  to  a,  b,  and  f  if 
the  rate  of  convergence  is  essentially  governed  by  the  approxima¬ 
tion  properties  of  the  trial  spaces  we  employ  and  the  smoothness 
of  the  solution. 

We  will  say  that  a  method  allows  for  postprocessing  if  it 
allows  a  significant  improvement  in  accuracy  to  be  obtained  by  a 
local  procedure  involving  only  a  single  element. 
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The  concept  of  a  generalized  finite  element  method  must,  of 
course,  include  practically  all  known  FEMs  satisfying  a),  espe¬ 
cially  the  standard  FEM  mixed  methods,  methods  using  different 
text  and  trial  functions,  and  methods  in  which  the  shape  of  the 
trial  functions  is  governed  by  the  differential  equations,  as 
e.g. ,  in  methods  for  the  convection-diffusion  equation  (see,  e.g., 

[  3 ]  and  references  therein) . 

Section  2  introduces  notation  and  preliminary  notions.  Sec¬ 
tion  3  introduces  and  discusses  the  Generalized  Finite  Element 
Method  (GFEM) .  Section  4  analyses  one  special  type  of  GFEM  which 
is  closely  related  to  finite  element  methods  based  on  L-splines 
(L  here  referring  to  a  differential  operator) .  This  method  is 
very  appealing  theoretically  and  is  extremely  robust,  but  in 
practice  has  various  drawbacks.  Nevertheless  it  provides  valuable 
insight  into  what  we  could  expect  at  "best."  Section  5  analyses 
a  modification  of  the  method  of  Section  4  which  better  satisfies 
requirements  of  implementation.  Section  6  analyses  a  further 
modification  which  is  still  very  robust  and  is  easily  implement  - 
able.  In  Section  8  we  show  it  is  identical  (in  the  sense  of 
leading  to  the  same  approximate  solution)  to  a  certain 
mixed  method.  In  addition,  the  method  is  shown  to  be  very  recep¬ 
tive  to  postprocessing.  Section  7  discusses  briefly  the  standard 
FEM  (as  a  GFEM)  and  shows  that  it  is  unsuited  for  problems  with 
rough  coefficients.  Also,  the  standard  FEM  is  shown  to  be  iden¬ 
tical  to  another  mixed  method.  Thus  it  is  seen  that  various  mix.  d 
methods  (or  versions  of  the  mixed  method)  have  substantially  dif¬ 
ferent  robustness  properties.  In  Section  9  we  further  elaborate  1  he 
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method  introduced  in  Section  6  and  show  that  it  "reacts"  very  well 
when  the  roughness  of  the  coefficients  is  reduced.  More  precisely, 
we  show  that  changing  from  measurable  coefficients  to  coefficients 
with  bounded  variation  improves  the  rate  of  convergence — as  we 
should  expect  from  a  robust  method.  Section  10  contains  some 
illustrative  computations  and  Section  11  presents  a  summary  of  the 
conclusions  of  the  paper. 

Some  of  the  results  of  this  paper  were  announced  in  [  4  ]  . 
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2 .  Preliminaries 

Throughout  this  paper  we  will  deal  with  one  dimensional 
boundary  value  problems.  Thus  we  consider  the  basic  model 
Dirichlet  problem 


(2.1) 


Lu  b  - (au1 ) '  +  bu  =  f , 
u ( 0 )  =  u(l)  =  0. 


0  <  x  <  1 


If  b  =  0  we  will  denote  the  differential  operator  by  instead 

of  L.  The  coefficients  a(x)  and  b(x)  are  assumed  to  be 
measureable  and  to  satisfy 


0  <  a  <  a(x)  <  B  < 


0  <  b(x)  <  B, 


0  <  x  <  1 


0  <  x  <  1 , 


Let  I  =  I  (tgjt^)  =  {x  :  tQ  <  x  <  t^,  where  0  5  tQ  <  <  1 . 

_  V 

1(0,1)  will  often  be  denoted  by  I.  By  W  (I)  with  k  =  0,l,... 

°k 

and  1  <  p  <  oo  we  denote  the  usual  Sobolov  spaces.  Wp  I)  denotes 

V  , 

the  subspace  of  Wp(I)  consisting  of  functions  which  vanish  to¬ 
gether  with  their  derivatives  of  orders  less  than  k  at  tQ  and 

Ir  t 

On  Wp(i)  we  have  the  usual  norms  and  semi-norms 


f'jJ-o 


(j)||P  xl/P 

Lp(I)  ' 


1  5  P  < 


k,p,  I 


max  II  u 
V.0- j-k 


L„, ( I )  ' 


k  ,  p ,  I 


Lp ( 1 } • 
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We  will  also  use  the  negatively  indexed  space  W  ^(i)  = 
($q)I))'  with  the  norm 


,  „  l/i  uvdxl 

I  U  ||  -i  T  =  SUP  — IT— t; - 

_1,P'1  v€#'(I)  HVlll,q,I 


where  ^  ^  =  1.  We  will  mostly  work  with  p 

r  Si 

If  p  =  2  we  write 


a)  and  p  =  2 . 


Hk(I)  *  Wk(I), 


|u||k,I  *  1  U  k  ,  2 , 1 ' 


and 


^  I  k  f  x  u  k  ,  2 , 1 # 


For  r  >  1,  Pr(I)  will  denote  the  space  of  polynomials  of  degree 
less  than  or  equal  to  r  on  I. 

The  solution  of  (2.1)  will  be  understood  in  the  usual  weak 
sense  in  ft'*’ ( I )  . 

In  our  analysis  we  will  use  various  a  priori  estimates  for 
the  Dirichlet  problem.  We  now  state  these  easily  proved  estimates. 

Lemma  2.1.  There  is  a  constant  0  <  C(a,8)  <  00  such  that 

(2.2)  Hull  1  -  5  C  (a  ,  8)  IlLull  0  r 

for  all  u  for  which  Lu  (  ^2 ^  and  u(0)  =  u(l)  =  0,  and 

<2.3)  llulll,-,r  5 
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for  all  u  for  which  Lu  e  ( I )  and  u(0)  =  u(l)  =  0. 

Lemma  2.2.  There  are  constants  0  <  C(a,R)  <  «  and 
0  <  C(a,g)  <  ~  such  that 

(2.4)  C(ct,3)  HLu||_1^  <  ||  u||lff  <  C  (a ,  B)  )|  Lu  |J_1  ^  — 

for  all  u  €  ^(1)  i  and 

(2.5)  C  (a ,  8)  If  Lu  If  _  x  f0Of  x  s  iiu,,l,o0,f  5  C(a,  8)  II  Lull  _1  ^  j 
for  all  u  €  W ■*"  (I)  . 
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3 .  Generalized  Finite  Element  Methods 

We  will  now  consider  a  class  of  finite  element  methods  for  the 
solution  of  (2.1).  Let  A  =  { 0  =  xQ  <  <...<  x  =  1} ,  where 

n  =  n(A)  is  a  positive  integer,  be  an  arbitrary  mesh  on  Y  and 
set  hj  =  Xj  -  xj_i  an^  Ij  =  l(xj_i»xj)  for  j  =  l,...,n,  and 
h  =  h ( A )  =  max h  .  .  With  (2.1)  we  associate  the  bilinear  forms 


j  J 

Bq(u,v) 

and 

B1(u,v) 

and  the  linear  functional 

F(v)  = 

The  exact  solution  u  of  (2.1)  is 

oi  — 

u  S  H  (I)  satisfying 

B  (u,  v)  =  BQ  (u,  v)  +  B-j^  (u, 

Given  a  finite  dimensional  space 
mation  u  6  S  to  u  by 

(3.1)  BQ(us,v)  +  B1(u, 

The  usual  choice  for  S  is 

S  =  Sr  =  { ^  €  H1 (T)  :  ^|T 

where  r  =  1,2,...  .  Then  u<,  =  u 


rl 

au ' v 'dx 

0 


!  buvdx 
1  0 


1 

fvdx. 

0 

then  characterized  as  the  unique 


)  =  F (v) ,  for  all  v  €  H^(I). 


c  ft'*' (I)  we  define  an  approxi- 


,v)  =  F(v),  for  all  v  €  S. 


€  Pr (Ij) ,  j  =  1, . . . ,n) 
is  the  standard  finite  element 


approximation  to  u. 
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Remark.  Sometimes  the  condition  S, 


°  1  — 

H  (i)  is  relaxed 


and  the  bilinear  forms  and  are  defined  as  sums  of  integrals 

over  the  subintervals  1^  .  In  this  paper  we  will  always  assume 

IT  o  X  _ 

SA  0  >  H  (I) .  However,  the  extension  of  the  ideas  developed  in  this 
paper  to  two  dimensional  boundary  value  problems  will  utilize  such 
a  generalization.  This  will  be  addressed  in  a  forthcoming  paper. 

The  finite  element  method  based  on  sf  n  uses  local  basis 
functions,  namely  those  functions  in  ^  that  are  characterized 

by  nodal  values  together  with  r  -  1  values  associated  with  each 


We  can  in  fact  write 


*  Sio®§A,0 


where 


=  {<f>  f  S^0  :  t(x.)  =  0,  j  =  0 ,  .  .  .  ,  n}  . 

It  is  clear  that  dim  Q  =  n  -  1  and  dim  =  (r-l)n.  Q  is 

called  the  space  of  internal  modes  or  a  bubble  space  and  is  charac¬ 
terized  in  various  ways,  e.g.,  by  function  values  at  r  -  1  in¬ 
terior  points  of  each  1^  or  by  values  of  the  first  r  -  1  moments 
on  each  I ^  .  This  is  especially  important  when  r  is  large.  In 
this  case  hiearchical  elements  are  of  importance;  for  more  on  these 
elements  see  [5],  e.g.  Local  basis  functions  lead  to  a  "collection" 
process  for  the  construction  of  the  sparse  global  stiffness  matrices 
from  the  local  stiffness  matrices  (which  are  associated  with  the  indi¬ 
vidual  1^).  The  load  vector  is  constructed  in  a  similar  way.  Wo  see 
immediately  that  this  essential  feature  of  the  standard  finite  element 
method  is  preserved  if  we  replace  the  finite  element  equation  (3.5)  by 
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(3.2)  Bq(A1ua,81v)  +  Bx  (a2UA 'B2V)  =  F<83v)'  for  a11  v  *  SA' 0  ' 


where  A^  , 

i  =  1,2 

,  and  8  .  ,  i 

l 

=  1,2,3,  are  one-to-one  linear 

mappings  of 

sl  0 

into  H1^) 

with  the  property  that  for 

'P  (  sl ,  o  ' 

Ml. 

and  8^1 

depends  only  on  xp  |  I  .  For  sim- 

3 

j 

3 

plicity,  we 

will  in 

addition  assume  that 

(3.3) 

(A )  (x_. 

)  =  ip  (Xj )  ,  i  =  1, 2 

and 

(3.4) 

<B.i|i)  (x. 

)  =  ^  (*j)  »  i  =  1, 2, 3,  j  =  1,  . 

(3.3)  can  be  weakened  while  (3.4)  is  essential.  Such  a  weakening 
is  important  in  the  two  dimensional  case.  (Cf.  the  previous  Remark.) 

We  note  that  the  standard  finite  element  method  corresponds  to  the 
choice  =  8i  =  E  =  the  identity  mapping.  We  will  refer  to  the  method 
defined  by  (3.6)  as  a  generalized  finite  element  method  (GFEM) .  The 
method  will  be  denoted  by  the  vector  L  =  (A^Aj,  B^,  of  mappings, 

and  the  corresponding  finite  element  approximation  by  ufi(U . 
eA(L)  =  u  -  ua(L)  will  denote  the  error. 

Remark .  The  existence  and  uniqueness  of  u^d)  as  defined 
by  (3.2)  cannot  be  established  without  further  assumptions.  This 
question  will  be  settled  for  each  of  the  specific  generalized  finite 
element  methods  discussed  in  the  paper. 

Practically  all  finite  element  methods,  including  methods  for 
solving  diffusion-convection  problems  and  reduced  integration 
methods,  can  be  put  in  the  form  (3.2).  The  form  (3.2)  obviously 
preserves  the  possibilities  of  the  collection  process  for  the  con¬ 
struction  of  the  stiffness  matrix  and  preserves  the 
sparsity  of  these  matrices,  so  that  one  of  the  main  ingredients  of 
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the  FEM  is  unchanged  (see  Section  1). 

A  number  of  questions  arise  in  connection  with  generalized 
finite  element  methods: 

i)  What  is  the  optimal  choice  for  the  mappings  ,  8^  , 

i.e.,  what  is  the  best  choice  for  the  local  stiffness 
matrices? 

ii)  How  do  implementational  requirements  influence  the 

choice?  How  much  effort  should  be  spent  in  the  construc¬ 
tion  of  the  A^  ,  B ^  as  compared  with  effort  spent  in 
mesh  refinement? 

iii)  Which  choices  lead  to  robustness  of  the  method? 

iv)  What  is  the  accuracy  of  the  generalized  finite  element 
method? 

These  questions  will  be  addressed  in  the  following  sections. 


13 


4.  The  L ^-Finite  Element  Method 

In  this  section  we  introduce  a  version  of  the  GFEM  which,  al¬ 
though  it  is  rather  impractical,  gives  valuable  insight  into  the 
best  that  could  be  expected  from  a  method  for  solving  (2.1)  when 
the  coefficients  are  rough. 

In  addition  to  the  usual  finite  element  space  Q  intro¬ 

duced  earlier  we  will  use  the  spaces 


S*"  =  {<|>  €  H1(I)  :  Ltp  j  _  6  Pr-2(  I.)  for  each  j} 

j  3 

and  q  =  fl  H3- (I)  .  (If  r  =  1  we  require  that  LtJj  =  0  on 

each  Ij  .  )  We  see  that  dim  Q  =  dimS^  =  nr  -  1.  We  now  dis- 

.  ~  r 

cuss  several  properties  of  the  spaces  . 

Lemma  4.1.  Given  a  function  u  €  H1 (T)  there  is  a  unique 
I^(a,b)  =  I^u  satisfying 


1 AU  <  SA 

(4.1)  ^  IAu(Xj)  =  u(Xj),  j  =  0,1, ...,n 

f  Q 

(u-IAu)(x-x.  ^)  dx=0,  1=0, . . .  ,r-2,]  =  l, . . .  ,n. 


IAu  €  ft1(I)  if  u  6  ft1*!) 


Proof .  The  existence  and  uniqueness  of  I^u  will  be  estab¬ 
lished  if  we  show  that  u  =  0  in  (4.1)  implies  I^u  =  0.  Thus 

-  r 

suppose  u  =  0.  Then  from  (4.1)  and  the  definition  of  we 


have 
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0  =  IAu  Ld^uldx 

1  j 

=  [a((I  .u)  'J2 +b(I.u)2]dx, 

■  I 

j 

from  which  we  get  I^u  =  0. 

We  will  refer  to  I&u  as  the  S^-interpolant  of  u. 

Lemma  4.2.  For  u  €  H1(I),  Ifiu  €  S^Q  is  characterized  by 

r1  f1  r 

(4.2)  (1  u) 'v'dx  =  u'v'dx,  for  all  v  €  SA  0 

0  A  J0  ' 

or  by 

(4.3)  B(Iau,v)  =  B (u, v) ,  for  all  v  f  SA,Q. 

Proof.  The  characterizations  (4.2)  and  (4.3)  can  be  seen  by 
writing 


(u-I  u)  ’v’dx  =  L[  (u-I.u)  (-v")dx  +  L(u-I  u)  (x.)Jv' (x.) 
0  JJ  j .  u  JaJ  J 


B(u-I.u,v)  =  l.  (u-I.u)Lvdx  +  L (u-I.u) (x.)Jav' (x.) , 

^  J  J  J  ^  J  '  '  J  J 

j 

where  Jav'(x^)  =  av' (xT)  -  av'  (x^)  .  We  note  that  Jav'(x.)  is 

well  defined  since  Lv |  €  P  (I„)  and  therefore 

x  x 

av' €  H  (1^)  for  l  =  l,...,n,  which  implies  av'  has  left 

and  right  limits  at  each  node  x.  . 

3 

Remark .  (4.3)  shows  that  IA u  is  the  Ritz  projection  of 

u  6  H1(I)  onto  the  subspace  SA  Q  with  respect  to  the  form  B. 


Lemma  4.3. 


For 


=  0  or  1 , 


k  >  0,  and  r  >  1,  we  have 
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(4.4) 


II u  -  1  auII  £  ,  Y  C  (a  ,  3 )  l|  Lull  ^  j 


and 


(4.5)  ||u  -  rAU|,£,oofi  5  C(a,B)hy||Lu||k>oo>I  , 

where  y  =  min  (k+2-£  , r+1-2. )  ,  and  C(ot,3)  is  independent  of  u  but 
depends  in  general  on  a,  3,  l,  k,  and  r. 


Proof .  We  first  prove  (4.4)  with  l  =  0.  Suppose 


Lw  =  - ( aw 1 )  '  +  bw  =  F  on  I 

w  ( 0 )  =  w(l)  =  0 


where  a  £  a(x)  £  6,  0  £  b(x)  £3  on  I,  and  let  v  be  the  Ritz 

projection  of  w  onto  the  subspace 


SQ  =  {0  ;  L0  €  Pr"2(I),  4>(0)  =  4.(1)  =  0 } 


with  respect  to  the  form  B(u,v)  = 


(au'v'+buv)dx,  i.e.,  suppose 


v  €  S, 


B(v,<|»)  =  B  (w,  <J>)  ,  for  cji  €  S 


We  write  a,  b  for  general  coefficients  in  order  to  distinguish 

from  a,  b  in  (2.1)  and  a,  b  introduced  below.  Now,  using  (2.2) 

(2.4),  and  the  fact  that  all  norms  on  the  finite  dimensional  space 
r  —  2  — 

P  (I)  are  equivalent  for  2  £  p  we  have 
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(4.6) 


II  Lv|| 


0  t  p ,  I 


Cl|Lvl|_1  j 
C  (a,  6)  II  vl!  1  J- 


S  C(a,B)  llw|llfY 


5  C(a,B)  ||F||0  y 
s  C(«,6)IIFII0jPjT  . 


Now  consider  u  -  I^u  on  I.  .  For  any  function  v(x), 

—  —  ^  _  x-x . 

x  f  I  .  ,  let  v (x)  =  v(x),  where  x  =  —  -  €  I.  By  a  standard 

j 

scaling  argument 


(4.7) 

where  I  u 


|u  -  !au»o,i.  5  Ch.1/2»u  -  T  ull0_Y 


.  ~  r  .  — 

is  the  S  -mterpolant  of  u  where 


Sr  =  { <f>  :  L4>  =  -(a4>* )  •  +  h2&|>  €  Pr“2(l)  }, 

i.e.,  I  u  is  the  element  of  S  having  the  same  values  as  u 

at  0  and  1  and  having  the  same  moments  of  orders  up  to  r  -  2 
as  u  (cf .  (4.1)). 

From  (2.2)  we  see  that 


(4.8)  ||u  -  I  u  H  o ,  Y  ~  c  («rS)  ||  L  (u—  I  u)ll0  j. 

Write 

Tu  =  (Tu)1  +  (Iu)  2 


where 


JL(Iu)1  =  0  on  I 

]^(Tu)1  (0)  =  u  ( 0)  ,  (Tu)  x  (1)  =  u  ( 1 )  . 
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It  is  easily  seen  that 

Cl[u  -  ( Iu)  ^  =  L  u  on  I 

Mu  -  (Tu)1l  (0)  =  [u  -  (Tu)^  (1)  =  0 

and  that  (T  u) 2  is  the  Ritz  projection  of  u  -  (Iu)1  onto 
SQ  with  respect  to  the  form  B(u,v)  =  / q  (au' v'+buv) dx,  where 

Sq  =  {<!>  :  L4,  e  Pr_2(I) ,  <)>(0)  -  <j)(D  =  0}. 

Hence  from  (4.6)  with  w  =  u  —  (I  u <  v  =  ( I  u)  2  *  and  p  =  2 
we  obtain 

(4.9)  ||L(Tu)||0fY  =  ||L(Tu)2ll0fI 

<  C  (ct,  B)  II  L  U||  0/Y  • 

Combining  (4.7),  (4.8),  and  (4.9)  we  have 

(4.10)  || u  -  IAu||0>I  5  C(a,6)hy2||L  u 

£  C(a,B)hV2||Lu  |lu_2fj 

since  y  >  2  for  £  =  0.  Now 


II  u 

XAu||0,Ij 

=  ||  (u— 4> )  - 

IA(u‘<)>)ll0,Ij 

for  any 

such 

that  L<J>  |  j 

■i 

€  pr_2(Ij)  . 

Combining  this  with  (4.10) 

we  see 

that 

J 

II U  - 

XAU^  0 , 1 j  ~ 

C(a,6)h  y2 

||Lu  -  L^IIm_2 fj  • 

Since 

L  ? 

is  arbitrary  in 

Pr_2( I)  we 

have 

k 
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I|U  '  TAu|lO,I.  -  C<«'B)h*/2  inf  IlLu  -  QIIm_2#y 

j  Q^pr“2 

~  C,c.8»h^2|LS|u.2>T. 

Then  a  further  scaling  argument  yields 

(4.11)  II  u  -  IAu|l0  T  <  C  ( a ,  B )  hlJ  |  Lu  |  2  . 

/  j  M  '  j 

(4.4)  with  £  =  0  follows  immediately  from  (4.11). 

(4.4)  with  £  =  1  is  obtained  by  a  slight  modification  of 
the  above  argument.  We  finally  note  that  (4.5)  follows  from  a 
similar  argument  based  on  (2.3)  and  (4.6)  with  p  =  °°. 

Lemma  4.4.  There  is  a  constant  C(a,B),  independent  of  u 
and  A ,  such  that 

(4.12)  H  IAU,I1,T  -  C  (a  ,  $)  Hull  x  ^ 

and 


(4.13) 


'  IAU^l/»f I  5 


C(af  B)  Hulllf00f f 


Proof .  (4.12)  follows  directly  from  the  coercivity  of  B. 

We  turn  to  the  proof  of  (4.13). 

Let  z  6  W^(I)  and  suppose  w  satisfies 


Lw  =  -(aw')'  +  bw  =  Q  €  Pr  2 (I)  on  I 

w ( 0)  =  z (0) ,  w(l)  =  z(l) 


1  .  rl 

wx-’dx  =  zx-’dx,  j  =  0,...,r-2 

0  0 


?  a(x)  -i  6,  0  b(x)  <  B.  Write  w  =  w^  +  w2  where 


where  0  <  n 
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w1  (1) 


z(l) 


and 


Lw2  =  Q 
w2  (0)  =  w2  (1)  =  0, 


We  also  write  z  =  z^  +  where 


7.1  €  H1(I)  ,  Zl(0)  =  z  ( 0)  ,  Z;L(1) 
g ( z1  # <() )  =  0,  for  0  €  ft1(I) 


=  z(l) 


and 


z2  €  ft1(I) 


ftl 


where  B(u,v)  = 


B(z 2,<p)  =  B(z,0),  for  0  €  H  (I), 
/^(au'v'+buvjdx.  B(z^,z2)  =  0  and  hence 


(4.14) 


I  z2  II 1 ~  C  (a ,  8)  II  z II  x  J  . 


Observing  that  w2  is  the  Ritz  projection  of  z2  onto  the  sub¬ 
space 

{<t>  :  L0  =  Pr~2(I)<t>  ,  0(0)  =  0(1)  =  0} 

we  see  that 


(4.15) 

We  now  show  that 


I W2 II 1 ,  l"  5  C(a'  II Z 2  II 1  I  ’ 


lWlll,oo,I  S 


1,“,I 


«JTT~ 


(4.16) 


C  (a,  6)  ||  z  || 
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First  we  show  that 


(4.17) 


1  Wlli ! ,«° ,1  '  c(«/3)  !l ZH o F«> ' 


If  b  =  0, 


w1  (x) 


-  Z  (1)  — 2  (0) 


+  Z  (0) 


and  the  result  is  immediate.  Now  suppose  b  ^  0.  Let  v  satisfy 


LQv  =  -(av*) '  =  0 

v(0)  =  z (0) ,  v (1)  =  z  (1)  . 


From  (4.17)  with  b  =  0  we  have 


(4.18) 


'  I,*,  i  5  c(a,B)  IIzIIq^j 


Letting  y  =  w^  -  v  we  see  that 


Ly  =  -bv 
y  ( 0)  =  y  (1)  =  0. 


Using  (2.3)  and  (4.18)  we  obtain 


11  y  '■  1 , 00 ,  "T  £  C(o,8)||vll0^j. 

(4.17)  follows  immediately  from  this  result  and  (4.18) 
Next  we  show  that 


(4.19) 


|W2!il,°°,I 


C(ra,  8)  II z I 


1,»,I  ' 


Using  (2.3),  (2.4),  and  the  equivalence  of  norms  on  the  finite 

3T-"  2  — 

dimensional  space  P  (I)  we  have 
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(4.20) 


2  1 ,°° ,  I 


C  (a  ,  g)  ||Q|| 
C  (a  ,  B  )  l|  Q II 


0,“>,I 


-1,1 

C(a,B)  IIw2H1>y 


Combining  (4.14),  (4.15),  and  (4.20)  we  obtain  (4.19).  (4.16) 
follows  from  (4.17)  and  (4.19). 

Now  consider  u  -  I^u  on  I_.  and  recall  the  notation  in¬ 
troduced  in  the  proof  of  Lemma  4.3.  By  a  scaling  argument  we  have 


(4.21) 


u-I.uL  T  £  ChT^|u  -  Tu| ,  T 
A  1 , 00 , 1  j  3  1  '!,«>,  I 


where  I  u  is  the  Sr-interpolant  of  u  where 


Sr  =  { :  L4>  =  -(a  <j>' )'  +  h2b<(>  €  Pr  2  (1)  } 


From  (4.16)  and  (4.21)  we  obtain 


(4.22) 

From  (4.22)  and  (4.5)  we  have 


5  c<a,e)h;1l|u|i1^rT . 


(4.23)  lu*IAul  lfCO/I_.  S  I  (u-Y)+Ia(u-y)  l1>00fl  +|y|  II-IaH^co,!. 

-1 


£  C(a,B)[h..  llu-Y  II  x  ^  oo ,  Y  +  lYlhj] 


for  any  constant  y.  Now  select  y  =  /J  udx.  Then 


and 


lu  -  vlll,~,T  *  c|u|l,»,l 


ItI  5  |u|0,-,T 


Thus  from  (4.23)  we  get 
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(4.2“)  <  CC.BHhj  lu^  +  h^ul  0,„,T3 

A  second  scaling  argument  applied  to  (4.24)  yields 


u  ^  AU  ^  1  / 00 , 1  j  "  C(a,p)  [  |  u  i  lf00/I^  +  hjlulo,«>,Ij 


C  (a,  B)  It u| 


j 


and  thus 


(4.25)  ^AU^1,°°,I.  “  C  (a  ,  S )  I!  u||  ^  ^  ^  ^  j  .  . 

(4.13)  follows  immediately  from  (4.25). 

The  main  goal  of  this  section  is  to  consider  the  L^-method, 

IT  ql  — 

where  1^  =  (A, A, A, A, A)  ,  the  mapping  A  :  q  -*  H  (I)  being 


defined  by  Au  =  I^(a,b)u  for  u  6  Q 


It  is  immediate  that 


A  is  one-to-one  and  onto  S.  n  ,  Au | T .  depends  only  on 

A ,  0  x-j 

u|  .  Au (x . )  =  u(x.),  and  that  A  1w  =  lj(l,0)w  for  w  (  Q  . 
3  3  3 

The  existence  and  uniqueness  of  u^(L^)  is  clear'  We  see  that 
u  €  n  implies  Au(x.)  =  0  and  that 

A  f  U 

ASi,0  '  ASlo*A§4,0- 

This  fact  plays  an  important  role  in  implementation. 

Lemma  4.5.  There  are  constants  C  and  C('x,(«),  independent 
of  A,  such  that 


(4.26) 


An!1  — 

AU  1,1 


C(a,  p)  !:ul 


1,1  ' 


(4.27) 


1 I 


C  ( a ,  6 )  I!  u ! 


,°°,I  '  f°r  U  *  SA,0  ' 


and 
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(4.28) 


!!  A  1w  H  - 

■L  r  ± 


C  Uw  H 


1,1  ' 


and 


(4.29) 


IA_1w 


,  Y  -  C||  wl|  ,  7  , 

1 1°°  ,  I  11  1  ,=»  ,  I 


for  w  €  S 


A,0 


Proof .  (4.26)  and  (4.27)  are  essentially  restatements  of 

(4.12)  and  (4.13),  respectively.  (4.28)  and  (4.29)  are  special 
cases  of  (4.26)  and  (4.27),  respectively,  corresponding  to  a  =  1 
and  b  =  0 . 


We  present  the  error  estimates  for  the  (^-method  in  the 
next  two  theorems.  Let  e^U^)  =  u  -  Au^l^). 

Theorem  4.1.  For  H  =  0,1  we  have 

(4.30)  ||  6  (I  )  i)  __  5  C(a,  8)hvl,lf  II  - 

1  i ,  i  k  > 1 

and 


(4.31) 


!e  .  ( I--, )  !|  - 

A  1  H ,  oo ,  i 


C  ( a ,  B)hu!lf  i 


k,  °°,  I 


where  y  =  min(k+2-£,  r+l-£) . 

Proof.  It  follows  immediately  from  the  definition  of  the 
(^-method  and  from  (4.3)  that  Au^i^)  =  IAu.  Thus  (4.30)  and 
(4.31)  follow  directly  from  Lemma  4.3. 
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Theorem  4.2.  For  £  =  0,1  we  have 

(4.32)  ,,eA  (Z-1}  (l£  ,T  s  C(a,6)h1_£|!f||0fI 
and 

(4.33)  l|eA(Ll)ll».,-fI  S  C(a#B)h1”a'Hfll0#<B>rT  . 

Proof.  These  results  follow  immediately  from  Theorem  4.1, 
the  fact  that  u^d^)  is  the  Q-interpolant  of  Au^d^),  a 

standard  approximation  result  for  Q  ,  and  (4.26)  and  (4.27). 

Remarks . 

1)  The  rates  of  convergence  given  in  (4.32)  and  (4.33)  are 
clearly  the  highest  that  could  be  proved  under  the  governing  hypo¬ 
theses  . 

2)  Suppose  that  we  have  found  u^d^).  Then,  although  the 
accuracy  of  uA  d^)  is  low  (see  Theorem  4.2  and  Remark  1), 

Aua(L1)  has  high  accuracy  (see  Theorem  4.1)  and  can  be  calculated 
locally,  element  by  element.  Thus  the  highly  accurate  Au,  (L.^) 
can  be  constructed  by  post  processing. 

3)  Theorem  4.1  shows  that  we  can  obtain  the  same  rate  of 

convergence  in  the  norms  II  •  !|  i  j  ,  £  =  0,1,  for  problems  with 

rough  coefficients  as  for  problems  with  smooth  coefficients  when 
post  processing  is  applied.  This  shows  the  maximum  possible 
robustness . 

4)  It  is  easy  to  show  that  postprocessing  does  not  improve 

the  rate  of  convergence  if  a  and  b  are  smooth  and  so  e,(L.) 

/ 1 
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and  e^U^)  are  of  the  same  order  in  this  case. 

5)  The  (.^-method  is  closely  related  to  Ritz  method  using 
L-splines  (see  [6],e.g.)  and,  as  would  thus  be  expected,  yields 
the  exact  solution  at  the  nodes. 

6)  We  have  presented  estimates  in  only  the  norms 

||  •  ||  ^  Y  anc^  II  *  II  £  j  f°r  £  =  0  and  1.  Parallel  results  hold  in 

the  norms  H*iljifp#i'  2sps°°,  £  =  0,1. 
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5.  The  L  -  Finite  Clement  Method 

III  this  s.i  H'  I  i  i  in  we  » ‘  t » 1 1 : :  i .  1 « •  l  I  1 1  <  •  1  •  ill<  •!  lit  til ,  wll<  ■  l  <  • 

L0  ~  ^A0'A0,A0'  A0,A0*  '  the  maPPin9  '^o  be;i-ny  constructed  from 

ic  r 

L  as  A  was  from  L  in  Section  4,  i.e.,  A^i  =  l^(a,0)u  for  u  <■  S . 
The  existence  and  uniqueness  of  u^d^)  is  immediate. 

Theorem  5.1.  For  2.  =  0  or  1  we  have 

(5.1)  "'’-VWU.I  5  Cla.eih^'llfll^j 

and 

(5.2)  I|u-*0UAU0)II*,”,T  "  C(a,S)h2-t|lfllki„iI. 

Proof.  Let  F  =  f  -  bu.  Then  from  (2.2)  and  (2.3) 

(5.3)  l,Fll0,I  S  C(ct'6)  l'f H0,T 

and 

(5.4)  »F»0,.,T  5  C<C,'6)I|£|I(,,«.J- 
Clearly 

(5.5)  (  Lqu  =  F 

\  u  (0)  =  u  (1)  =  0. 

Let  u,  (L  )  denote  the  -  finite  element  approximation  to  the 
solution  of  (5.5).  From  Theorem  4.1  applied  to  (5.5),  (5.3)  and 

(5.4)  we  have 

(5.5)  =-A0S4U0M,j  5  C(.'>,B)h2',ilf:  0_J 
and 


*  < 
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(5.7) 


Writing 


(5.8) 


5  C(a,B)h2'Vl!0i„jT. 


W  =  W  +  z(t0> 


we  immediately  see  that 

(5.9)  B ( AgZ , Aqv)  =  B(A0UA(L0) ,A0V> 

-  VVi'V’V*  -VVi'V’V* 

'  BllU-A0“4<L0>'A0V)'  £°r  V  f  S4,0- 

From  (5.9)  and  the  coercivity  of  B  we  have 


il A o“ I! i , i  "  c 65 B (A0z,Aoz) 


<  c(.,6)iiu-A0uAa0)iig ;Tiia0,ii0iT 


and  this,  together  with  (5.6),  yields 


(5.10) 


I! Aoa  11 1 , Y  5  c(a'B)h  llfll^j- 


(5.1)  follows  from  (5.6),  (5.8),  and  (5.10). 

For  l  =  0,  (5.2)  follows  from  (5.7),  (5.10)  and  the  fact 

that 


(5.11) 


^A0Z^0,°°,I  5  l|A0z||l,  I* 


It  remains  to  prove  (5.2)  for  l  =  1.  It  is  easily  shown 
that  there  is  a  positive  constant  C(a,S)  such  that 


! 
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(5.12) 


|Bq (u,v) | 

SUP  irrm - “  -  C(a,6)||v|L  .  T>  for  u  €  w;(I) 

Ol  —  llnlooT  ±  t  J.  ,  ±  1 

u€WA(I)  '  ' 


From  (4.3),  (4.13),  and  (5.12)  we  have 


(5.13) 


iup 


Bq (w, v) 

fwT 


,  r  " "  11 1  co  T 
W€A0SA  ,0 


=  sup 

ufwrT) 


|  Bq  (IA  (a,  0)  u,v) 
i|lA(a,0)u|l1^fI 


=  sup 
U?W* (I) 


Bq(u,v) 

lu|li,«.T 

'1,1,1' 


U 


1,°°,! 


ZA(a'0)ul'l,-,T 


>  C(a,0)  ||v|L  .  T,  for  v  €  A.S 


r 
0“A 


(5.13),  together  with  the  fact  that  an  operator  and  its  adjoint  have 
equal  norms, yields 


(5.14)  sup 


Bq(w,v 


VfA0SA,0 


- —  >  C(a,8)||w||  j, 

1,1,1  ,X 


for  w  €  A0Sa. 


Now  from  (5.7),  (5.9),  (5.10),  (5.11),  and  (5.14)  we  obtain 


(5.15)  II A  z ||  m  -  <  C  (a,  8)  sup 


0  "!,«>, I 


V  €S 


A ,  0 


=  C(a,8)  sup 
,r 


LvviVj 


lBl<U~V  U0)-A0z,40v) 


ves 


TV*!,!,! 


A ,  0 


“  C(a,6)  ( II  u-A0ua  (  Lq  )  H0fOOfI  +  HAozl'o,"o,  i} 


s  C  (a,  8)  n  ||f  || 


0,oo,I 


(5.2)  with  Z  =  1  follows  directly  from  (5.7),  (5.8),  and  (5.15). 
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Remarks . 

1.  The  rate  2—2.  in  (5.1)  and  (5.2)  cannot  be  improved. 
This  can  be  seen  by  examining  the  special  case  a  =  1  and  b 
rough . 

2.  (5.2)  with  Jl  =  0  shows  that 

|  u(Xj  ) -uA  (Lq)  (Xj)  I  5  C  (ct ,  6 )  h2  II  f  II Q  j  =  l,...,n-l. 

We  also  have  the  analogue  of  Theorem  4.2. 

Theorem  5.2.  For  1  =  0  or  1  we  have 

(5.16)  l|eAU0)ll^,I  5  C(a,B)h1_Jl||f||0fT 
and 

(5.17)  llW't,-,?  5  C(a,8)h1-t||f!l0i„>T. 

Remark .  Comparing  Theorems  5.1  and  5.2  we  see  that  the 
Lq  and  methods  have  the  same  rates  of  convergence.  On  the 

other  hand, implementation  of  the  LQ  method  is  much  simpler 
since  the  equation  LQ  ip  =  <j>  can  be  solved  by  quadrature. 
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6 .  The  L 2  -Finite  Element  Method 

In  this  section  we  consider  the  =  (A g,E,E,E,E)  method, 

where  E  is  the  identity  mapping  and  is  as  defined  in 

Section  5.  First  we  analyze  the  case  b  =  0. 

Lemma  6.1.  Suppose  b  =  0.  Then  for  l  =  0,1  we  have 


(6.1) 

and 

(6.2) 

where  u 


5  C(a,8)h1J||f!lk>T 

lu-Vi11!1 11  >.,«■, T  5  C(o.,6)hu||f!|k_<0_1 

min (k+1 , r) . 


2 

Proof.  We  will  prove  (6.1)  in  detail  for  r  =  2.  In  S.  n 
-  A  9  u 

we  construct  a  basis  by  choosing 


satisfying 


,4> 


n-1 


n 


W  *  5ij-  i,j  ■ 1 . .  - 1> 

f1  4>.dx  =  0,  i  =  l,...,n  -  1,  j  =  1 

j  1 


n. 


C^(Xj)  =  0,  i  =  l,...,n,  j  =  l,...,n  -  1, 

/  j  C^dx  —  5  ^  ’  3  =  1 »  •  •  •  r  n . 

“1  -J 


For  f  €  Hk(I)  let 
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1 

1 


where  6  is  the  Dirac  distribution  at  x.  and  X.  is  the 

X  3  3  3 

characteristic  function  of  I..  It  is  easily  shown  that 

3 

A0ua(12)  is  the  exact  solution  of  problem  (2.1)  with  right 
hand  side  fA.  We  note  that  this  fact  implies  the  existence  and 


uniqueness  of  u  (Lj) 


Thus  from  (2.4)  we  have 


(6.3) 


'u  A0UA  (  L2)  "l,  I  S  c(a,  6) 


Now  let  w  be  the  solution  of 


(6.4) 


(  -w"  =  f,  x  i  I 
|w(0)  =  w(l)  =  0 


and  let  w^  be  the  standard  finite  element  approximation  to  w 

2 

in  the  subspace  SA  Q  i.e.,  let  wA  be  the  Ritz  projection  of 

2  1 
onto  S  n  with  respect  to  the  form  /_  u'v'dx.  As  above,  we 

f  U  U 

see  that  wA  is  the  exact  solution  of  (6.4)  with  right  hand 
side  f  .  Thus  from  (2.4),  applied  to  the  problem  (6.4),  we 
obtain 


(6.5)  5 

Combining  (6.3),  (6.5),  a  well  known  estimate  for  the  standard 

finite  element  method,  and  a  simple  regularity  estimate  for 
(6.4)  we  have 
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liu-A,0UA  ^  L2*  "l,  I  ~  C  (a ,  0)  h  llwHk+2>I 

S  C(a,0)h2||f||k 

This  is  (6.1)  for  r  =  2.  The  proof  for  general  r  is  similar. 
Finally  we  note  that  the  proof  of  (6.2)  follows  the  same  lines 
but  is  based  on  (2.5)  instead  of  (2.4). 


Remark .  In  general  one  cannot  obtain  a  better  rate  of 
convergence  in  H^(I)  (W^(I))  than  one  obtains  in 

H1(I)  (wJ(T)). 

Lemma  6.2.  Suppose  b  =  0.  Then  for  l  =  0  or  1  we  have 

(6.6)  l|eA  (  L2}  !l  Jl  5  C(a,0)h1_eif||o>1 

and 

(6.1)  l|eA(L2)|li,I  5  C(a,0)h1"^||f||OfOOfI. 

Proof .  The  proof  of  this  lemma  is  similar  to  the  proof  of 
Theorem  4.2. 

We  now  turn  to  the  general  case. 


Theorem  6.1.  There  is  a  constant  hQ  =  hg(u,8)  such  that 
for  0  ■-  h  <  hg ,  u^dg)  is  uniquely  determined.  Furthermore 
there  is  a  constant  C(a,8)  such  that 


(6.8) 

and 

(6.9) 


leA(L2)!l0fT 


C(U,  B)hllf  llg  Y 


|eA  (  L2)  ^  0 , 00 , 1 


C  (a,  0)  h||  f  i 


0 ,  I 


O’ 


for  0  <  h  <  h 
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Proof.  Consider  the  problem 


-(aw' ) '  =  F,  x  f  I 
w{0)  =  w(l)  =  0. 


Let  T  be  the  corresponding  solution  operator,  i.e.,  let 
TF  =  w,  and  let  TA  be  the  L 2  approximate  solution  operator, 
i.e.,  let  TaF  =  wA(I-2).  We  regard  T  and  Th  as  operators  on 
L2(T).  From  Lemma  6.2  we  have 


(6.10) 


lT-T4liL2,L2  5  C(a,B)h. 


Let  u  be  the  exact  solution  of  (2.1)  and  u  (L  )  the  i-2 


approximation  to  u.  Then 


(6.11) 


(6.12) 


u  =  T  (  f  -  bu) 


UAU2) 


Vf  "b  UAU2M 


Using  (6.11)  and  (6.12)  we  obtain 


(6.13)  u  -  uAU2)  =  (T-TA)f  -  Tb(u-uA(L2)) 

+  (Ta-T) b(uA ( L2) -u)  +  (TA-T)bu 


which  can  be  written  as 


(6.14) 


[  I  +  Tb  +  (T  -T)b]eA  U2)  =  (T-r[V  (f_bu) 


where  I  is  the  identity  operator  on  L2(I).  Now  TbI  is  com¬ 
pact  on  L2  and  it  is  easily  seen  that  I  +  TbI  is  one-to-one. 
Thus  (I  +  Tbl)  ^  exists  and  is  bounded.  From  (6.10)  we  see  that 
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l|  (Ta-T)  bl  ||  L  L  <  C(a,B)h 

and  hence 

(6.15)  l|  [I  +  TbI  +  (T-T)bI]_1|L  T  s  C (a , 3) 

A  L2'L2 

for  h  sufficiently  small.  Identity  (6.14)  together  with  the 
invertibility  of  [I  +  Tbl  +  (T^-Tjbl]  implies  the  existence  and 
uniqueness  of  uA(l0).  It  follows  immediately  from  (6.10),  (6.14), 
and  (6.15)  that 


!|eA  (  L2}  ,!0 , 1  "  C(a»B)h||f-bu||0  y 

S  C  (a ,  3)  h||  f  ||  0  ^  — . 

This  is  (6.8).  The  proof  of  (6.9)  is  similar. 

Remark. 

1.  The  assumption  that  h  is  small  is  not  necessary  as  will 
be  shown  in  Section  8.  See  Remark  2  following  Theorem  8.2. 

2.  In  Lemma  6.2  the  case  l  =  1  is  essentially  a  stability 
result,  this  result  is  also  valid  in  the  general  case  (when  b  ^  0) . 
See  Remark  2  following  Theorem  8.2. 

Theorem  6.2.  For  i  =  0  or  1  we  have 
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(6.18)  I!u-A0uA  a2)  h,I  -  C(a,B)hllf|l0  j 
and 

(6.19)  Hu-A0ua(L2)  11  lt~ti  -  C(af  6) h|| f  II o ^ oo # Y" 

Proof .  Recall  the  notations  from  the  proof  of  Theorem 
6.1.  Now  let  be  defined  by  T^F  =  A^T^F.  Just  as  we 

obtained  (6.13)  we  now  obtain 

(6.20)  u  -  A0ua(L2)  =  (T-TA)f  -  Tb(u-uA(L2)) 

+  (TA-T)b(uA(L2)-u) 

+  (TA-T)bu. 

From  Lemma  6.1  we  have 

(6.21)  llT-T.D  .  S  C(o,B)h||f||n 

A  L2,H  0,1 

Using  (2.2)  we  see  that 

(6.22)  ||Tb(u-uA(l2)  )  !|  y  5  C(a,  B)  ||u-uA(L2)  ||Q 

Finally,  combining  (6.20)— (6.22)  and  Theorem  6.1  we  see  that 

llu  -  V  <L2>»1,T  5  C«»,S)hllf|l0iT. 

This  gives  (6.18).  The  proof  of  (6.19)  is  similar. 

Remarks . 

1.  Estimates  (6.6)  and  (6.8)  cannot  be  improved  when  general 
a  and  b  are  present.  If  a  and  b  are  smooth  then,  of 
course,  better  estimates  can  be  derived,  namely 
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(6.23)  ileA<L2)  !I£,I  S  C!|f||k,YhU- 

where  <  =  0  or  1  and  u  -  min(k+2-v',  r+l-O.  Here  C 
depends  on  the  smoothness  of  a  and  b  as  well  as  on  a  and  8. 

Similar  results  hold  in  connection  with  the  norm  ||  •  |]  — . 

£  ,  00  ,  I 

2.  (6.16)  and  (6.17)  canno'  in  general  be  improved. 

Theorem  6.2  shows,  however,  that  the  accuracy  can  be  improved 
by  post  processing. 

3.  The  rate  of  convergence  in  (6.23)  is  the  best  possible, 
as  can  be  seen  from  the  theory  of  n-widths  (see  [  7  ] ) .  if 

a  and  b  are  merely  measurable  it  is  also  possible  to  show 
that  (6.8)  and  (6.16)  give  the  best  possible  rate  of  convergence. 
Thus  we  see  that  u^(l2)  yields  the  optimal  rate  of  convergence 
in  all  cases.  Similar  results  hold  in  connection  with  the 
norm 

4.  The  12  method  is  well  suited  for  implementation.  In 

fact  we  will  see  in  Section  8  that  it  can  be  obtained  as 

a  mixed  method  or  can  be  easily  implemented  directly. 

5.  We  could  also  discuss  other  choices  for  the  and 

8^  which  would  modify  in  various  ways  the  accuracy  of  the 
approximate  solution. 
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7.  The  Standard  Finite  Element  Method  (The  L^-method) 

The  standard  finite  element  (or  Ritz)  method  corresponds  to 
the  choice  1 ^  =  (£»£,£,£,£)•  In  this  section  we  discuss  a 
negative  result  that  shows  the  contrast  with  the  ^-method.  We 
assume  r  =  1. 

Consider  the  family  of  boundary  value  problems  given  by 

a  ,  b  =  0,  f  =  1,  where  A  =  {x_,...,x  }  =  { 0 1  ' 
n  n  u  n  n  n 

and 


r 


2.  Xj.j 


X  .  , +x . 

<  X  <  -1-i— 3 


an(x)  =  < 


1, 


X  .  ,  +x  . 

3-1  .3  < 

2 


x  <  x..  ,  j  =  1 ,  .  .  .  ,n 


and  let  un  denote  the  exact  solution.  It  is  easy  to  see  that 

u^  ( L 3 )  (the  standard  finite  element  approximation  to  u  )  is 
n 

the  piecewise  linear  interpolant  to  the  exact  solution  of  the 
problem 


From  Lemma  6.2  we  have 


n 


-  u. 


n 


a2) 


0,1 


c 

n 


and  it  can  be  shown  (see  Theorem  9.3)  that  u^  converges  in 
L2  to  the  exact  solution  of 
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f-3w"  -  1 

I  w(0)  =  w(l)  =  0. 

This  shows  that  the  L  approximation,  (L 2 ),  is  very 

n 

accurate  while  the  standard  finite  element  approximation, 

uA  (L^),  fails  to  "converge",  i.e.,  un  -  — *■  0. 

n  n 

Remarks . 

1.  We  see  that  the  standard  method  is  very  nonrobust.  It 
does  not  converge  for  rough  coefficients  at  all,  and  gives  very 
poor  results  when  a  changes  significantly  between  the  nodel 
points.  In  contrast  the  L method  performs  well  whether  the 
coefficients  are  smooth  or  rough.  Because  the  implementation  of 
the  /_2  involves  the  same  amount  of  work  as  does  the  standard 
method,  the  ^-method  should  be  prefered  in  all  cases. 

2.  We  have  analyzed  only  linear  boundary  value  problems. 

Many  of  the  ideas  introduced  here  are  valid  for  nonlinear  problems 
also.  For  strongly  nonlinear  problems  the  L 2  method  is  very 
promising . 


1  nr  '  rn 
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8 .  The  Relation  Between  Generalized  Finite  Element  Methods  and 
Mixed  Methods 

In  this  section  we  discuss  the  relation  between  Generalized 
Finite  Element  Methods  and  mixed  methods.  Mixed  methods  for  the 
approximate  solution  of  (2.1)  can  be  derived  from  the  displacement 
formulation  (2.1)  by  introducing  an  auxiliary  variable,  writing 
the  second  order  scalar  equation  as  a  first  order  system,  casting 
the  systems  into  variational  form,  and  then  discretizing  the 

resulting  variational  equations.  This  can  be  done  in  various 
ways.  We  consider  two  mixed  methods  for  (2.1). 

Letting  s  =  pu 1  we  write  (2.1)  as  a  first  order  system 
in  two  ways,  namely  as 
/■ 

au'-s  =  0,  0  <  x  <  1 

<  -s'  +  bu  =  f,  0<x<l 

^  u ( 0 )  =  u(l)  =  0 

and 

r  S 

u' - =  0,  0  <  x 

a 

\  -  s'  +  bu  =  f,  0  <  x 

u  ( 0)  =  u(l)  =  0, 


<  1 

<  1 


and  then  consider  the  associated  variational  formulations 


T 


40 


r 


u  f  H1 (I) ,  a  €  H° (I) 


(8.1) 


au'adx  - 


0 


sodx  =  0,  for  all  o  €  H°(I) 


1  fl 

sv'dx  +  I  buvdx 
0  -1 0 


fvdx ,  for  all  v  €  ft1 (I) 


and 


r  U  i  H1 (I)  ,  '  f  H° (I) 


(8.2)  < 


fl 

^  so 

u'adx  - 

—  dx  = 

J 

1  0 

0  3 

r1 

r  1 

sv'dx  + 

buvdx  = 

V  i 

1 0  J 

0 

fvdx,  for  all  v  €  H1(I). 


0 


The  mixed  methods  we  consider  are  now  obtained  by  discretizing 
(8.1)  and  (8.2)  : 


(8.4) 


r 


UA  f  SA  ,  0 '  SA  ( 


(8.3)  < 


au^odx  - 


s^v’dx  + 


s.odx  =  0,  for  all  a  €  d 


bg  vdx 


ua  6  sa.o'sa  ( nr1 


fi 

o 


<  J„UAadx  - 


sa 


s^v'dx  + 


bu^vdx 


fvdx,  for  all  v  e  S 


A  ,0' 


—  dx,  for  all  a  e  S 


r-1 


fvdx,  for  all  v  €  S*"  n, 
0  A  '° 
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where  S^-1  =  f ip  €  H°(I):  ip/].  =  Pr-1(I..),  j  =  l,...,n). 

As  indicated  above,  the  mixed  formulations  (8.1)  and  (8.2) 
are  obtained  from  the  displacement  formulation  (2.1)  via  the 
introduction  of  an  auxiliary  variable.  Clearly  this  process 
can  be  reversed,  i.e.,  (2.1)  can  be  obtained  from  (8.1)  or  (8.2) 

by  eliminating  the  variable  s.  It  is  also  possible  to  eliminate 
this  variable  at  the  discrete  level,  i.e.,  to  eliminate  s^ 
from  (8.3)  or  (8.4),  obtaining  a  Generalized  Finite  Element 
Method.  This  is  made  precise  in  the  next  two  theorems. 


Theorem  8.1.  The  elimination  of  s^  from  (8.3)  leads  to 

the  L -  finite  element  method  for  u. ,  i.e.,  u.  is  the 

3  A  A 

standard  finite  element  approximation  to  u. 


Proof.  Suppose  u^ ,  satisfies  (8.3).  Letting  a  =  v' 

in  the  first  equation  in  (8.3)  we  find  Jgau^v'dx  =  J^s^v'dx. 
Using  this  in  the  second  equation  we  see  that 


r 


UA  €  SA,0 


au!v'dx  + 
,  A 


bu . vdx 
,  A 


[  fvdx,  for  all  v  €  sf 

Jo  A' 


i.e.,  u^  is  the  L finite  element  approximation  to  u. 


Theorem  8.2.  The  elimination  of  from  (8.4)  leads  to 

the  L_  finite  element  method  for  u. . 

2  A 

Proof.  Suppose  uA,  s A  satisfies  (8.4).  From  the  first 

equation  in  (8.4)  we  see  that 
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u 


i  i 

AHj 


=  P 


I  . 
3 


r  “  1 

where  P  is  the  L„-projection  onto  P  (I.)  .  The  mapping 


3A|l.  "  Rj(SA  I.1 
1  3  J 


=  P 


-1 


i.i  from  pYi1.) 


to-one;  let  Qj  =  •  Then 


r-1 

to  P(I  )  is  one 


sa 1 1 j  *  Vily 


Using  this  in  the  second  equation  in  (8.4)  we  see  that 
satisfies 


UA  €  SA  ,  0 


yn 

zj=i 


_  Q.u!  v'dx  +  bu.vdx 
I.  ]  A  \n  A 


fvdx,  for  all  v  (  S 


0 


It  remains  to  show  that  Qju^  j 


=  a(A0UA) '  I. 


From  (4.2) 


see  that  P(AQ<|>)  '  |  T  =  «p '  | :  for  «p  €  S^Q  and  from  the 


~  y  XT  ”  l 

definition  of  we  have  (A^u)'^  *  P  ( I j ) -  Thus 


V^WIl-j3  *  P<*oVilJ 


J 


=  u 


a|i3 


we 


This  completes  the  proof. 


M  <! 


Remarks . 


1.  If  r  =  1  both  the  L ^  and  L methods  can  be  viewed 
as  arising  from  the  application  of  the  standard  finite  element 
method  to  a  problem  with  appropriately  perturbed  coefficients. 
Consider  the  problem. 


-  (a^u ' ) '  +  bu  =  f ,  0  (  I 


u (0)  =  u(l)  =  0 


where  a^  is  constant  on  each  1^,  and  consider  its  approximation 
by  the  standard  finite  element  method.  If  a^  =  a^ ,  where 

^ I . adx 

a^ I  =  — ^ -  ,  i.e.,  if  a^  is  the  piecewise  average  of  a, 

Tj  j 

we  obtain  the  standard  finite  element  method  for  the  solution 
of  (2.1).  Also,  it  is  clear  (when  r  =  1)  that 


Q 


I . 
3 


I . 


’  dx 

I  .  a 
3 


where  Q ^  is  the  operator  introduced  in  the  proof  of  Theorem  8.2, 


Hence  if  a.  =  a. ,  where  a. , _ 

A  A  A  I  . 


r  dx  V 
'  I  .  a 

_J _ 

h. 

3 


-1 


,  i.e.,  if 


is  the  piecewise  harmonic  average  of  a,  we  get  the  method. 
It  is  obvious  that  and  a^  can  differ  significantly  if  a 

changes  significantly  on  the  subintervals  I.,  while  for  smooth 

—  2 
a,  a^  and  a^  differ  by  0(h  ). 

The  observation  that  the  L 2  method  is  tha  same  as  the 
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standard  method  applied  to  the  problem  with  a  replaced  by 
a  ,  where  a^  is  the  piecewise  harmonic  average  of  a,  shows 
that  the  L ^  method  can  be  easily  implemented,  its  implementation 
being  of  the  same  difficulty  as  for  the  standard  method. 


2)  Using  the  operator  R^  introduced  in  the  proof  of 

Theorem  8.2  we  can  prove  the  existence  and  uniqueness  of  uA (L^) 

for  all  h.  See  Remark  1  following  Theorem  6.1.  R.  is  positive 
r— . 

definite  on  F  (I.)  with  respect  to  the  L_-inner  product  and 

3  £ 

its  eigenvalues  are  bounded  away  from  00  and  0  uniformly  in  j. 


Thus 


f  2  r-1 

(Q.z)zdx  t  X  z  dx,  for  z  6  Pr  A(I.) 

JI.  ]  I  .  3 

3  3 


with  X  >  0  independent  of  j .  Now  suppose 


C1  r 

at AgU^ ( L 2 ) ] 1 v '  +  buA(i2)v  dx  =  0,  for  all  v  €  SA  Q. 

Then 

o  -  EQju'(L2»u'(L2)  d*  ♦  |Wu2)J2 

j 

-  AIj|  tuA(L2)]2dx 

Tj 

which  implies  uA(i2)  =  0* 

3)  Nemat-Nasser ,  in  [8]  and  a  subsequent  series  of  papers, 
suggested  the  use  of  the  variational  principle  (8.2)  in 
conjunction  with  trigonometric  polynomial  approximating  functions 
for  the  approximate  calculation  of  eigenvalues  of  problems  with 
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rough  coefficients.  The  authors  derived  error  estimates  for 
this  method  and  corresponding  finite  element  methods  (cf.  (8.4)) 
in  [9]. 
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9 .  Additional  Properties  of  the  L ^ -method 

So  far  in  this  paper  we  have  assumed  the  coefficients  a(x) 
and  b(x)  are  in  L^d).  In  this  section  we  derive  additional 
properties  of  the  L  ^-method  in  the  case  in  which  a(x)  has  bounded 
variation  and  r  =  1.  We  begin  by  proving  a  perturbation  of  coef¬ 
ficient  result. 

Let  u  denote  the  solution  of  (2.1),  i.e.,  let  u  be  the 
a  a 

solution  of 


/'-(au^)  1  +  bua  =  f,x  €  I 
lu  (0)  =  u  (1)  =  0 

a.  3. 


and  let  u  be  the  solution  of  the  corresponding  problem  with  a 
aA 

replaced  by  afiI  where  is  the  piecewise  harmonic  average  of 

a  subordinate  to  the  mesh  A  =  { 0  =  <  x2  <  •  •  • 

the  step  function  defined  by 


xn  =  1),  i.e. 


Hi. 


dx 


All. 


h. 

3 


We  now  present  an  estimate  for  u  -  u  under  the  hypothesis 

a  a^ 

that  a(x)  has  bounded  variation . 

For  a  function  c(x)  with  bounded  variations  denote  by  cA 
the  piecewise  average  of  c.  We  will  assume  functions  with  bounded 
variation  to  be  right  continuous  at  each  x  <  1  and  left  contin¬ 
uous  at  1.  It  is  easily  seen  that  V(cA)  <  V(c),  where  V  =  V-- 
denotes  the  variation  on  I. 


Lemma  9.1.  There  is  a  constant  c  such  that 
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(9.2) 


x  . 

f  3  A  . 

(c-c  )  <Pdx 


0 


2-1 


5  CV-(c)  I  4>l ,  _  ^  h 


l.p.I 


for  <P  €  W  (T)  ,  1  S  p  £  00 ,  j  = 


Proof.  Let  4>  €  W^(I)  »  let  be  the  piecewise  average 


of  <t>,  and  set  $(x)  =  /q  )dt.  Then 


(9.3) 


x  . 

f  3  A  A 

(c-c  )  <Pdx 


x  . 
f  3 


A  ,  ,  A 

(c-c  )  ( <P-<P  )  dx 


x  . 
3 


(c-c  ) $ 1 (x) dx 


(c-c  ) d $ 


X  . 

3 


<J>d  (c-cA) 1 


5  CIUIL  —  V(c-CA) 
0 ,  I 


ClUlt  -  V  (c)  . 
0  ,«>,  I 


Now  we  bound  ||$|L  — .  Suppose  x .  .  s  x  <  x  .  Then 

U,0°,I  Jfr—  1  X. 


*(x) 


li-l 


(<p -<P^)  dt  +  j  (4>-4>A)dt 
Jxa-i 
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fX 


( <J>-4>  )  dt 


Vl 


and  thus  by  the  Poincare"  inequality  we  have 


(9.4) 


P-1 

»(x)|  5  »tIU,|'o,p,It  ht  P 


(9.2)  follows  from  (9.3)  and  (9.4). 

For  b  =  0  we  obtain  two  estimates  for  u  -  u 

a  aA 

is  an  estimate  for  ||u  -a,,  ||n  _  =■  and  the  second  is  an 

a  , p  t 1 

for  u  (x.)  -  u  (x.),  j  =  1 , .  .  .  ,n-l . 
a  J  3a  J 


Theorem  9.2.  Suppose  b  =  0.  There  is  a  constant 
depending  only  on  a  and  3  satisfying 


(9.5) 


( u  — u  )  zdx  <  C  (a  ,  S )  V—  (a )  ||  z | 
I  a  aA 


^  '  a  aA  ' - 1  *  -  I  ,!“"0,q,Il|fl|0,q,I 


for  z  €  Lp(I)  ,  1  <  p  5  00,  where  ^  +  i  =  1. 


The  first 
estimate 

C  (  a  ,  8  ) 

1  + 
h 


Proof.  The  solution  ufl  of  (9.1)  is  given  by 


r-g  a, 
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where 


■  I 


X  /q  f ( s) ds 


a  (t) 


dt 


and 


E(x)  = 


dt 

a  (t)  ' 


and  the  solution  is  given  by 


(9.7) 


u  (x)  =  -D  -(x)  + 

aA  aA '  * 


D  -(1)E  .  , 

aA'f  aA(x) 

(1) 

aA 


We  first  estimate  [n(D  £-D  ,)zdx. 

J  0  a ,  f  a  .  .  r 

A 

Set  Q (x)  =  /*  z(t)dt.  Then 


r 


(D  f-D  f)zdx  = 
Q  a , r  aA,t 


:( 


fl/fX  x  2.  rt 

jQ  (HTtr~aA(t))j0 


f  (s)dsdtjz(x)  dx 


■  1\  A 


Noting  that  —  =  (— )  we  can  apply  Lemma  9.1  with  c  =  — , 

aA  a  a 

<t>(x)  =  /q  fdsQ(x),  and  x^  =  1  to  get 


(9.8) 


(D. 


_f-D  )zdx  5  CV,i)|*|1;q;Th1+P  Ul0  f 


1+- 


5  CM  VU)llfll0  q  f  llzll0>qrTh  P. 
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Next  we  estimate  /g  (E-E^)zdx.  We  see  that 


(E  -E  )zdx  = 
0  a  aA 


fX  1  1 

( -  -  —  )dt  z  (x)  dx 
0^0  3  3  A 


0  ^WQ<’t,d!'' 


Now  we  apply  Lemma  9.1  with  c  =  —  and  <Mx)  =  Q(x)  and  obtain 

3 


(9.9) 


(E  -E  )z  dx  £  C(«)VT(a)  I1  z  t!  - 

g  a  aA  J ' q ' x 


l+i 

h  P 


1  r  X 

Applying  Lemma  9.1  with  c  =  --  and  *(x)=/Q  f(s)ds  we  see  that 


(9.10) 


D  (1)-D  f (1) !  - 
a , r  a&  '  r 


1  1  1  fX 

(i-_L)  f  (s)dsdt  | 
0  3  A  ■'0 


£  C(a)vT(a)  Hfll0  qrT  h 


From  (9.6)  and  (9.7)  and  the  fact  that  E  (1)  =  E  (1)  we 


see  that 


(9.11)  u  (x )  -  u  (x) 


D  (x)  -D  Jx)  +  E  (1)  [D  (1)E  (x)  -D  f(l)E  ,  . 

aAff  a'f  3  3'f  3  3A'£  3A(x) 

(D  _(x)  -D  ,(x)]  +  E(1)_1{[D  ,(1)  -D  -  (1)]E  (x) 

aA,f  a , f  a  a,f  a 


+  D  -(1)  [E  (x)  - E  (x) ]  ) 
a A,r  a  3 A 


a*.’  r  ^ 
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(9.5)  follows  immediately  from  (9.11)  and  (9 . 8) - (9 . 10) . 


Theorem  9.2.  Suppose  q  t  0.  Then  there  is  a  constant 


C(a)  such  that 


(9.12) 


ua(Xj)  -u  <x.)  I 


1+- 

C(a)VT(a)  l|fl|0#q#T  h  P 


Proof.  From  (9.11)  we  see  that 


u  (x.)  -  u  (x. 


a  3 


aA  3 


-1 


D  _(x.)  -  D  ,(x.)  +  E  (1)  [D  f(l)  -  D  f ( 1 ) ] E  ( X . ) 
a^,r  ]  a,r  j  a  a,r  aA,z  a  ^ 


X  x 

Applying  Lemma  9.1  with  c  =  -  and  st(x)  =  f  f(s)ds  we  obtain 

d  u 

(9.12)  . 


We  now  consider  the  general  case  in  which  b  /  0, 


Theorem  9.3.  There  is  a  constant  C(a,P)  such  that 


1+i 


(9.13) 


u  -  u 


7  s  C  (a »  6)  V  -(a)||f!jn  „  =■  h 


a  'a.  "  0  ,  p,  I  -  . i'-'"‘"0,q,I 


for  1  5  p  5  oo. 


Proof.  If  b  =  0  we  can  apply  (9.5)  to  get 


|ua  "  UaAll0,p,I 


sup  _  • 
z*L  (I) 

A 


^  ( u  -u  )z  dx  | 

U  3  3  . 

A 


0  E  I 

'p-1' 


C(a)  V  (a)  ||  f 


l+± 

.  P 


0  -P  T 
U'p-1'1 
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If  we  let 


LqU  =  - (au' ) ' , 


LA,0U  -  -<aAU,)'- 


and 


Bu  =  b(x)u 


this  result  can  be  written  as 


(9.14) 


-1  -1  1+o 

lL0  “LA,0^l  (I)  ,L  (I)  5  C(ot)  V(a)h  , 

^3  P 


where  —  +  —  =  1,  and  the  desired  result  amounts  to  proving 
P  ^3 


l"L0+B)_1  -  aA,0+B)'i|lLq(T),Lp(T)  5  C(„)  V(a)h 


1  + 


We  observe  that 


(9.15)  (L.+B)  1  -  (L+B)  1 


(LA+B)  1(L+B) (L+B)"1  -  (L+B)"1 


(La+B) _1LaLa1 [ I-B (L+B)  h  -  [-(LA+B)  1B+I]L  1L(L+B)  1 
[I-(LA+B)"1B]LA1[I-B(L+Br1]  -  [- (La+B) _1B+I ] L_1 [I-B (L+B) _1 ] 


[I-(LA+B)_1B]  (L^-L-1)  [I-B  ( L+B)  _1  ]  . 


T)  I  M 
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It  is  easily  seen  that 


(9.16) 


-1. 


i-o.4+b>  bii  (T)>  ,r)  <  cfa.d 


and 


(9.17) 


I-B(LtB)-\  j  (I  <  ClcBI 
q  q 


It  now  follows  immediately  from  (9. 14) - (9. 17)  that 


|ua-“aA»0,p,T 


(L4+b)_1<-  (I,+B>'l£»0,p,T 


II  -  (La+B)-1B|IL  _L  l|L-1-L-1HL  (T)iL  (I)  III  -BIL+B)"1!^  >L 
P  P  ^  P  ^  ^3 


0  ,  q ,  I 


1+- 


5  C  (a  ,  0)  Vy  (a)  II  f  II 0  ( q  f  Y  h 
This  completes  the  proof. 


We  now  return  to  a  consideration  of  the  L^-method  under  the 
assumption  r  =  1. 


Theorem  9.4.  There  is  a  constant  C(ot,0)  such  that 


(9.!8)  „U  "  uAl,0rPrr  -  C(ct,  6 )  v^(p)||fj|0  rq,T  h 


1+— 
P 


Proof.  Using  Theorem  9.3  we  have 
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(9.19) 


U  'UA  O,p,I  *  +  l|UaA'UAl!0,p,I 


l+i 


s  C(a,3)VT(a)h  P  llfll0,q,i  +  Hua  “Vo^.T 


We  next  observe  that  u.  can  be  viewed  as  the  standard  finite 

A 

element  approximation  to  u  (cf.  Remark  1  following  Theorem  8.2). 

a  A 

If  b  =  0,  is  the  piecewise  linear  interpolant  of  u  and  by 

A 

standard  approximation  results  we  have 


(9.20) 


0,p,I  5 


1+- 


C  (a ,  6)  h  p  ||  f  I 


0,q,I 


If  b^o,  (9.20)  can  be  proved  by  a  slight  modification  of  the 
usual  duality  argument.  (9.18)  follows  directly  from  (9.19)  and 
(9.20)  . 


Remark .  The  estimate  given  in  (9.18)  is  optimal. 

Theorem  9,5.  Suppose  b  5  0.  Then  there  is  a  constant  C 
such  that 


(9.21) 


u(x.)  -  uA(Xj)  | 


C(a)Vj(a) 


f|l0,q,T  h 


1+- 
P 


for  1  <  p  5  °°. 

Proof.  This  result  follows  immediately  from  Theorem  9.2  and 

the  fact  that  u.(x.)  =  u  (x). 

A  ]  aA '  D 
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Remarks . 

1)  Theorem  9.4  gives  the  rate  of  convergence  of  the 
i-2-method  without  post  processing.  It  is  possible  to  prove  that 
the  same  result  holds  for  the  L^-method. 

2)  One  can  prove  that  under  the  hypotheses  of  Theorem  9.4 

we  have  the  rate  of  convergence  h1//p  in  the  norm  IML  __  y  • 

•L  $  P  •  x 


; 


i 


i 

i 

i 

i 


.  t 
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10 .  Some  Illustrious  Computations 

In  this  section  we  illustrate  the  improvement  ot  ten-il  t.y  i  he 
L2-method  over  the  L^-method  tthe  standard  1 liute  element  method)  by 
considering  two  examples. 


Example  1.  Consider  problem  (2.1)  with 


a  (x) 


a^  ,  0  <  x  <  1/2 

ar  ,  1/2  <  x  <  1 


b  (x)  =  0 


where  a^  and  ar  are  positive  constants  and  consider  the  approxi¬ 
mation  of  the  solution  u  by  the  (^-approximation  u^  (L2)  and 

n 

the  (.^-approximation  ( L corresponding  to  the  uniform  mesh 

n 

with  an  odd  number  of  subintervals  (x^  =  jh,  j  =  0,1,..., n, 
h  =  n  ,  n  odd) .  One  can  obtain  exact  expressions  for  the  errors 
at  the  nodes.  For  the  node  just  to  the  left  of  I^x(n-l)/2^ 
we  find 


eA  (L3)  ( 2  “  2} 
n 


1 1  h  i  / 1  y  /  X  h » 

u(2  ”  2]  ~  UA  (L3}  2  ~  2^ 
n 


h(ar~a£} 
8a£ (a£+ar) 


eA  (  L 2 }  { 2  ~  2 ) 
n 


h  (a£-ar} 
8a£(a£+ar) 


2  .  2, 


2  2, 


{8afar(a‘+ap  +  (a^-ap 


u  (2  “  UA  (  L2^  (2  “  2} 

n 


h  <arar} 

8a£(a£+ar) 


h  (ar-a£) 
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From  these  expressions  we  see  that  if  ar  -  a?  is  large 
relative  to  h,  then  e^iL^)  <<  e(L^),  while  if  a^  -  af  is 

_  u  ^ 

approximately  equal  to  h,  then  eA(L2 )  eA^L3^  - j.  Tables 

16a* 

1  and  2  give  numerical  values  of  the  errors  at  the  nodal  points 
for  the  case  a^  =  1,  ar  =  100. 


X 

eA  <L3> 
n 

.2 

9.23-3 

-1.96-3 

.4 

1.85-2 

-3.92-3 

.6 

-1.85-4 

3.92-5 

.8 

-9.23-5 

1.96-5 

Table  1. 

eA  (L3)  and  e A  ( L2 ) 
n  n 

ar  =  100,  n  =  5. 

for  a£  = 

X 

eA  U3) 
n 

eA  (L2> 
n 

1/21 

5.32-4 

-2.65-5 

2/21 

1.06-3 

-5.29-5 

3/21 

1.60-3 

-7.94-5 

4/21 

2.13-3 

-1.06.4 

5/21 

2.66-3 

-1.32-4 

6/21 

3.19-3 

-1.59-4 

7/21 

3.72-3 

-1.85-4 

8/21 

4.25-3 

-2.12-4 

9/21 

4.79-3 

-2.38-4 

10/21 

5.32-3 

-2.65-4 

11/21 

-5.31-5 

2.71-6 

12/21 

-4.78-5 

2.44-6 

13/21 

-4.25-5 

2.18-6 

14/21 

-3.71-5 

1.92-6 

15/21 

-3.18-5 

1.65-6 

16/21 

-2.65-5 

1.38-6 

17/21 

-2.12-5 

1.11-6 

18/21 

-1.59-5 

8.35-7 

19/21 

-1.06-5 

5.58-7 

20/21 

-5.30-6 

2.84-7 

Table  2.  eA  U3) 
n 

and  eA  (L  ) 
n 

for  a£  =  : 

ar  =  100,  n  =  21. 


T 
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Example  2 .  Here  we  consider  an  example  with  a  smooth  but 
rapidly  varying  coefficient.  Let 

500 ( x  -  h 

a(x)  =  6  +  - y~,  b  (x)  =  0, 

1+100|  x-~| 

and  consider  the  same  mesh  as  in  Example  1.  To  give  some  idea  of 
the  variation  of  a\x)  we  present  a  table  of  its  values  (Table  3). 


X 

a  (x) 

.0 

1.10 

.1 

1.12 

.2 

1.16 

.3 

1.24 

.4 

1.45 

.5 

6.00 

.6 

10.55 

.7 

10.76 

.8 

10.84 

.9 

10.88 

1.0 

10.90 

Table  3.  Values  of  a(x). 


Table  4  gives  the  errors  at  the  node  just  to  the  left  of  1/2 
for  the  standard  method  (L3)  and  the  L2 -method  for  this  example. 


n 

e4  U3> 
n 

eA  {L2) 
n 

9 

1.48-3 

-4.31-4 

15 

6.31-4 

-1.15-4 

21 

3.34-4 

-4.60-5 

27 

2.00-4 

-2.28-5 

33 

i . 30-4 

-1.28-5 

39 

9.03-5 

-7.75-6 

? 

45 

6.53-5 

-5.12-6 

4. 

eA  {L3}  and  SA  U2] 

at  the  node  just 

- 

n  n 


to  the  left  of  1/2. 
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11 .  Conclusions 

In  this  section  we  summarize  in  tabular  form  some  of  the  re¬ 
sults  we  have  obtained  on  the  convergence  rates  for  the  various 
methods  and  then  make  several  remarks  on  conclusions  that  can  be 
drawn.  The  convergence  rates  in  the  table  are  for  smooth  f. 


-method 

(standard 

FE  method) 

1^ -method 

L 2 -method 

lleV 

Hell  , 

H 

_ 

Hell  Q 

H 

llell  , 

H 

Hell  Q 

H 

Hell  . 
H 

a,b  smooth 

hr+1 

■ 

hr+1 

m 

hr 

a,b  measurable, 
without  post¬ 
processing 

h° 

h° 

h 

h° 

h 

u0 

h 

a,b  measurable, 
with  post¬ 
processing 

h° 

hr+1 

hr 

h 

h 

a  bounded 
variation, 
b  measurable,  r=l 

h° 

h° 

h3/2 

h3/2 

h1/2 

Remarks . 

1)  The  standard  FEM  is  very  nonrobust. 

2)  The  rates  for  the  L^-method  show  that  by  investing  a 
considerably  larger  effort  in  the  construction  of  the  micro¬ 
stiffness  matrices,  one  can  significantly  increase  the  accuracy- 
-as  compared  with  the  standard  method.  This  suggests  that  the 
goal  of  minimizing  the  entire  computational  effort  with  respect 
to  desired  accuracy  could  well  be  served  by  increasing  the  per¬ 
centage  of  total  effort  that  is  spent  at  the  micro  level.  This 


could  be  especially  important  in  convention  with  adaptive  pro 
cedures . 

3)  The  computational  effort  involved  in  implementing  the 
L^-method  is  exactly  the  same  as  with  the  standard  method.  T 
the  l2~method  should  clearly  be  preferred  over  the  standard 
method  in  all  cases. 

4)  The  Generalized  Finite  Element  Method  offers  a  much 
larger  freedom  in  the  choice  of  computational  procedures  than 
does  the  standard  finite  element  method,  and  as  a  onsequence 
offers  the  possibility  of  significant  improvement  ±n  accuracy 
especially  when  used  in  conjunction  with  adaptive  procedures. 
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